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In dynamical spacetimes apparent and event horizons do not coincide. In this paper we propose a geometrical 
measure of the distance between those horizons and investigate it for the case of the Vaidya spacetime. We show 
that it is well-defined, physically meaningful, and has the expected behaviour in the near-equilibrium limit. We 
consider its implications for our understanding of black hole physics. 
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I. INTRODUCTION 

If you pull on a rubber band, it stretches. As the Moon or- 
bits Earth, its gravitational field raises the ocean tides. Grav- 
itational waves impinging on an interferometer signal their 
presence by causing its arms to stretch and shrink. In each 
of these examples a physical object is distoited by external 
forces, and that distortion signals a transfer of energy. Intu- 
itively, one would expect black holes to behave in a similar 
way: tidal forces and gravitational waves should distort them 
and a distorted black hole should throw off gravitational waves 
as it rings down to equilibrium. Indeed numerical simulations 
and perturbative calculations show such behaviours (see, for 
example, |1 1 and [2] respectively). However at a fundamental 
level things are not so simple. 

First, the physical objects discussed in the opening lines 
are locally defined: one can easily identify the edge of a rub- 
ber band or the surface of the Earth's ocean. This contrasts 
strongly with a classically defined black hole. A causal black 
hole is a region of spacetime from which no (causal) signal 
can escape. Its boundary is the event horizon. For outside 
observers, this null surface is the boundary between the un- 
observable events inside the black hole and those outside that 
can be seen. While very intuitive, a little thought shows that 
this definition is necessarily teleological: one determines the 
extent (or existence) of a black hole by tracing all causal paths 
"until the end of time" and then retroactively identifying any 
black hole region (FIG.[T]i. Defined this way a black hole is a 
feature of the causal structure of the full spacetime rather than 
an independent object in its own right. 

The second complication follows from the first: a black 
hole does not directly interact with its environment. By defini- 
tion it is causally disconnected from the rest of the spacetime 
and while it can be affected by outside events, signals from its 
interior cannot escape to influence the surrounding spacetime. 
Again this contrasts quite strongly with our initial examples: 
one can see and feel a stretched rubber band and signals from 
the Earth can certainly reach the Moon. 

The third complication also follows from the teleological 
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FIG. 1: A schematic showing the various horizons for a spherically 
symmetric spacetime with the angular dimensions suppressed. Hori- 
zontal location measures the areal radius of the associated spherical 
shell while the direction of increasing time is roughly vertical outside 
the event horizon but tipping horizontal-and-to-the-left inside. The 
shaded gray region is a shell of infalling null dust, gray dashed lines 
represent "outgoing" radial null rays, and the apparent horizon is that 
associated with a spherically symmetric foliation of the spacetime. 



definition of black holes: event horizons evolve in unexpected 
ways. As an example, consider FIG.fTlagain which depicts a 
simple spherically symmetric black hole which transitions be- 
tween two equilibrium states by absorbing a series of shells of 
infalling matter Intuitively one might expect the event hori- 
zon to remain unchanging until the first matter arrives, grow 
while it falls in, and then return to equilibrium when all of the 
matter has been absorbed. This is not what happens: the event 
horizon is expanding before the matter arrives and its advent 
actually curtails, rather than causes, further expansion. 

In spite of these apparent problems, it has already been 
noted that in practical calculations event horizons do seem 
to reflect physical interactions between a black hole and its 
surroundings in the way that one might expect. What then is 
happening in those examples? Well, as is recognized by the 
membrane paradigm |3|, it is the external near-horizon grav- 
itational fields ("just" outside the event horizon) that are the 
agents of interaction with the environment. From this per- 



spective the importance of horizons is not that they directly 
interact with the rest of the world, but rather that they might 
closely reflect the physics of those near-horizon fields. For ex- 
ample, in an interaction with gravitational waves, a shearing 
of a horizon should indicate a similar shearing of the gravita- 
tional fields "just-outside" the horizon. Further, a radial inflow 
of matter into a black hole will increase the strength of those 
fields and result in a corresponding increase in horizon area. 
More topically, recent work has linked horizon deformations 
to "kicks" in black hole mergers [4J. 

Though this all seems plausible, we have already noted that, 
in general, event horizon evolution does not reflect local in- 
fluences in such an obvious way. Understanding the regime 
where they will evolve in the expected way is one goal of this 
paper In pursuit of this, we recall an alternate way to define 
black holes which is geometric and based on trapped surfaces: 
closed and spacelike two-surfaces with the property that all 
normal congruences of forward-in-time travelling null rays 
decrease in area. These surfaces are characteristic of black 
hole interiors. For example, a classic theorem by Penrose|l5| 
tell us that if the null energy condition holds, then those nor- 
mal congruences necessarily terminate at spacetime singular- 
ities. Further in asymptotically flat spacetimes, trapped sur- 
faces are necessarily contained in causal black holes|6, 7| . 

Given a foliation of spacetime into "instants", one can (in 
principle) locate all trapped surfaces at each instant and take 
their union to find the trapped region. The boundary of 
that union is an apparent horizon and can be shown to be 
marginally trapped with vanishing outward null expansion|6 
|7). Now in practical calculations, such as in numerical rel- 
ativity, one cannot possibly locate all trapped surfaces. In- 
stead one works with the outermost marginally trapped sur- 
face (which can be identified by well-known algorithms LSJ). 
Even without an underlying spacetime foliation, one can study 
hypersurfaces foliated by marginally trapped tubes as poten- 
tial black hole boundaries. In general these are marginally 
trapped tubes but specific types of them have been much 
studied over the last couple of decades including trappingf^l, 
isolated|10|, and dynamical horizons lfTTI . 

Now, these geometrically defined horizons can be locally 
identified and, just as importantly, they evolve in response to 
local stimuli. For example, in FIG. [T] the geometric horizon 
expands if and only if matter is falling through it. This makes 
them attractive as candidates for understanding black hole dy- 
namics. However, there are still difficulties: 1) in asymptot- 
ically flat spacetimes on which the energy conditions hold, 
marginally trapped surfaces are necessarily contained within 
event horizons and so still not in causal contact with the rest of 
spacetime and 2) marginally trapped tubes are non-rigid and 
may be deformed |12-14|. That is, unlike event horizons, they 
are non-unique. 

So from a theoretical perspective, both the causal and ge- 
ometric definitions are problematic if we hope to understand 
black holes as astrophysical objects interacting with their sur- 
roundings in intuitive ways. However, whatever theory says, 
in perturbative calculations, one sees both event [2| and ge- 
ometric 1 15 1 horizons interacting with their surrounding in 
the way one might expect: shearing and expanding under 



the influence of gravitational waves. Perturbative calcula- 
tions necessarily describe near-equilibrium black holes, and 
in this regime the various theoretical problems somehow melt 
away. Indeed in the perturbative calculations, apparent and 
event horizons even coincide to a high degree of accuracy 1 2 1 . 

We then advocate the following model of near-equilibrium 
black holes. By definition, the near horizon fields of the mem- 
brane paradigm are just outside the event horizon which in 
near-equilibrium should also be very "close" to the event hori- 
zon. In this regime we would expect that both the apparent and 
event horizons would closely reflect the evolution of the near- 
horizon fields and so be a useful tool for understanding black 
hole evolution. The goal of this paper is to (begin to) quantita- 
tively understand the details of how this happens. In particular 
we will develop tools to geometrically quantify what "close" 
means. Though our ultimate goal is to develop these ideas in 
general, in this paper we will focus on spherically symmetric 
spacetimes and in particular the Vaidya solution as an initial 
testing ground. Some similar ideas have recently been advo- 
cated and developed in lfT6l and will be discussed at appropri- 
ate points of this paper 

We proceed as follows. In Section II we review the mathe- 
matics of trapped surfaces and horizons, specializing to spher- 
ical symmetry which is all that we will need for this paper 
In particular we consider the circumstances under which ei- 
ther event or apparent horizons may be considered to be near- 
equilibrium. Section III turns to Vaidya spacetimes, exactly 
locating their apparent horizons. In the appropriate regimes 
event horizons are located either numerically or perturba- 
tively. In both cases we again pay particular attention to the 
near-equilibrium regime and note that under those conditions 
the horizons are close as expected (at least in a coordinate 
sense). Section IV quantifies that observation by introducing 
a physically motivated, geometrically defined notion of the 
(proper time) distance between the horizons. We explore and 
refine this with numerical and perturbative calculations and 
prove that for near-equillibrium Vaidya spacetimes, the event 
and apparent horizons are indeed arbitrarily close together. 
Section V reviews the preceding sections and examines their 
further implications. 



II. MATHEMATICS OF SPHERICALLY SYMMETRIC 
HYPERSURFACES 

As promised, we now review the mathematics of horizons: 
that is the geometry of two- and three-dimensional surfaces in 
a four-dimensional spacetime. Since this paper focuses on the 
Vaidya spacetime, we restrict our attention to spherically sym- 
metric spacetimes and (hyper)surfaces that share those sym- 
metries. For details of derivations or the corresponding non- 
symmetric expressions see, for example, lfT3l . 



A. Geometry of two-surfaces and three-surfaces 

Let {M,gab) be a time-orientable four-dimensional 
spherically-symmetric spacetime. Then given a spacelike 



two-sphere S (which shares the spacetime symmetries), the 
normal space at each point has signature (1 + 1) and can 
always be spanned by two future-oriented null vectors. We 
respectively label the future-oriented outward and inward 
pointing null normals as £ and n, cross-normalize so that 
£ ■ n = — 1 and require that they also share the spacetime 
symmetry: if ^ is any rotational Killing vector field then 
^ r£ ^ ^ rn ~ 0. There is a degree of freedom for these 
vectors: they will remain null, cross-normalized and with 
the correct symmetries under rescalings i —> h£ and n -^ n/h 
where h is any function that shares the symmetries of S. 

The induced metric on S is that of a sphere and may be 
written in standard spherical coordinates as 



qAB == r^{[d9]A ® [de]B + sin^ 0[<i0]A ® [d<j)]B) 



(1) 



where r is the areal radius. The corresponding area element is 



^yq — r sin , 
and the (two-dimensional) Ricci scalar is 
- 2 



(2) 



(3) 



Next, consider the extrinsic geometry. Then, the symme- 
try again significantly simplifies matters, and the only non- 
vanishing components of the extrinsic curvature are the ex- 
pansions: 



0(0- 



Qi,i\ — 



^/q r 



and 



(4) 



Two-surfaces can be classified by the signs of these null 
expansions. If n is inward pointing we expect 0(„\ < (areal 
radius gets smaller moving inwards) and classify the surfaces 
as: 

1. untrapped if din > 0, 

2. marginally trapped if 9/(\ — and 

3. trapped if difn < 0. 

As mentioned in the introduction, the interiors of stationary 
black holes are made up of trapped surfaces while their bound- 
aries may be foliated with marginally trapped surfaces. 

Now consider a hypersurface H foliated by spacelike two- 
spheres. Then for some function C and an appropriate assign- 
ing of the null vectors, we may always write a tangent vector 
to H as 



r = r-cn''. 



(5) 



The sign of C determines the signature of the H: C < <=> 
timelike, C = 4=> null and C > <^ spacelike. The rate of 
change of the area element "up" H is 



^ 



rVl' 



q [di^ii) ■ 



C0, 



\n)) 



(6) 



We will refer to C as the evolution parameter. 

There are several other derivatives "up" the surface that are 
also of interest. First we have 



Ky = -n„y"Vj' 



(7) 



As suggested by the notation, this reduces to the surface grav- 
ity (or equivalently the inaffinity of the scaling of the null vec- 
tors) for a Killing horizon. Next the derivatives of the expan- 
sions are 



■^r%) =^rd{i) - l9^f^/2 + Gah£"£ j 



(8) 



^a_,Jj 



C ( ^ - Gabi"n ^ - 9(() 0(„ 



and 



^ r 0(„) = - Kr e^„)+C e' II + Gabrfn 



,M,J) 



(9) 



+ 



fia_,J} 



+ Gabtn — 0(^) 0(„) 



These follow directly from Equations (2.23) and (2.24) of 
IIT3I (taking into account our spherical symmetry which sets 
shears, normal connections, and surface derivatives to zero). 
The 1/r^ comes from a two-dimensional Ricci scalar term. 



B. Spherically symmetric event horizons 

An event horizon is a null surface generated by null 
geodesies. In terms of the mathematics of the last section, 
it is a three-surface with C = so that 



if^%)=jc^%)- 0(-^^,)/2 + G, 



lah 



tt 



(10) 



This is the Raychaudhuri equation and it encapsulates the 
counter-intuitive behaviours of event horizons. Given the Ein- 
stein equations and null energy condition and adopting an 
affine scaling of the null vectors (iCf = 0), it follows that 



■^iQ(() < 






(11) 



The rate of expansion of a null surface naturally decreases 
with time. There is nothing profound about this statement. 
The expansion as defined by (pi is defined relative to the size 
of the sphere at the time of calculation. Thus the rate of ex- 
pansion of a sphere expanding with constant speed naturally 
decreases, even in flat space. However, a non-zero flux of 
positive-energy density matter will cause the rate to decrease 
even more. Again this is to be expected: the gravitational in- 
fluence of more mass inside the shell should decrease the rate 
of expansion. 

The second law also follows directly from Eq. (111. By 



this equation if Oi(\ is ever negative for a congruence of null 
curves, then that congruence will necessarily converge to a 
caustic. However, no event horizon can contain a caustic |6|. 
Thus O/n > and event horizons are non-decreasing in area. 



For a stationary event horizon Ou-j — and so at first thought 
one might try to define a near-equilbrium horizon to be one 
with 6u\ <C l/r. However this is, of course, not a geometri- 
cally meaningful statement since the null vectors retain a seal- 



Further, since J^yOin = 0, C can be found from dSl as 



ing freedom. Instead we return to ( 10 1 and motivated by the 
thermodynamic analogy define: 

Definition: A spherically symmetric null surface is said to 
be slowly evolving if 



Jf f;0(^) < K(.e^e) and -Q^,, « G„i,r^^ 



(12) 



Note that this is a scaling-invariant statement: a quick check 
will show that if we rescale £ ^ hi and n — ^ n/h, then the net 
effect is to rescale both conditions by a factor h^ which will 
not affect the inequalities. 

If these conditions hold then, as first discussed in 1171, it 



is conventional to reinterpret Eq. (10 1 as a first law of event 
horizon mechanics with the (apparently) counter-intuitive be- 
haviours only showing up at higher order. We find 



Ke^(^= ^qKiQ(e) ~ V^Gahiy^ , 



(13) 



which is of the same form as the T5S — 5E near-equilibrium 
version of the first law of thermodvnamics OOl . Achieving 
this form is the main motivation for our definition. It will 
be further supported by later calculations but for now note 
that this first law is causal (to lowest order) with matter fluxes 
through the horizon driving the expansion. 

Next we consider geometrically defined horizons. 



C. Spherically symmetric geometric horizons 

A marginally trapped tube (MTT) is a three-surface that can 
be foliated with marginally trapped spacelike two-surfaces. 
There are several special cases of MTTs. Given the null en- 
ergy condition, a null MTT is a non-expanding horizon and 
physically corresponds to a black hole in equilibrium with its 
surroundings lITOl . Dynamical black holes are represented by 
dynamical horizons which are spacelike and expanding fTTl. 
Future outer trapping horizons (FOTHs) can be either non- 
expanding or dynamical and, in addition to 9f„\ < 0, have 
^ „9i(\ < (there are fully trapped surfaces just inside the 
horizons) [9|. Here we will review some of their basic prop- 
erties, but for a more general discussion see aforementioned 
original sources or one of the review articles lfT9l . 

With 0(£) = 0, the expansion ^ of the MTT is: 

■^■rV^ = -V^Ce^„) (14) 

and so with 0(„) < 0, 

1. C < 4^ timelike and contracting (^/y/q < 0) 

2. C = <=> null and constant area (J£ f^/q = 0) 

3. C > <=> spacelike and expanding (J£ y^/q > 0) . 



C = 



Gabi"e' 



\/r^--Gabi"nb 
For a marginally trapped surface 



■^(7%) ~ ^ ( ^ ^ GabPn 



(15) 



(16) 



and so for a FOTH the denominator of Eq. (15 i is positive 
by the assumption that there are fully trapped surfaces just 
inside the horizon. Violations of this condition only occur un- 
der extreme conditions such as the formation of new horizons 
outside of old ones Ii20l422l . 

Thus by (fTSJl we see that, if the null energy condition holds 
and Gi,bi"i^^ 0, a FOTH is necessarily a dynamical hori- 
zon (and vice versa). On the other hand if GabC-"^^ — then 
C — Q and the horizon is null, non-expanding and (weakly) 
isolated. The causal nature of these geometric horizons is 
clear: they expand when matter falls in and are otherwise qui- 
escent. Thus, like event horizons, they cannot decrease in area 
and so obey a second law. 

Finally we consider the conditions under which a geomet- 
ric horizon may be considered to be near-equilibrium. The 
conditions are more complicated than for an event horizon but 
are well-studied lfT3l [TSl |23| and simplify greatly in spheri- 
cal symmetry. Again we are guided by the requirement that a 
near-equilibrium horizon should obey a first law of the form 
( 13 1. To that end we combine (IHl) and 
obtain 



with Om = to 



fcr0(r) -C^r0(„) = Gah^i^-C^ (G„fc«V + 02,)/2)(.17) 

Then we define a slowly evolving MTT as follows. 

Definition: A spherically symmetric marginally trapped 
tube is said to be slowly evolving if 



if r0(„) « Kr0(„)| and C(0(;,)/2 + G„i«"«'') « ^ . (18) 

Again this is invariant under rescalings of the null vectors 
and with the help of ( 15 i is sufficient to reduce ( 17 i to a 755 w 
5E form of the first law: 



K^yB(r)~GabtP 



(19) 



This first law is not the only motivation for this set of con- 
ditions. For example 



1 



,ceL « 



1 



^,^r<l, 



(20) 



where y" is the unit-normalized version of Y". That is, the 
rate of change of area radius is small relative to distance mea- 
sured "up" the horizon. Further discussion and motivation can 
be found in fT3l . 

The non-uniquess of geometric horizons is not apparent 
from the above discussion. This is due to our spherical 
symmetry restrictions. For spherically symmetric spacetimes 



there is a unique spherically symmetric marginally trapped 
tube. The other marginally trapped tubes don't share those 
symmetries. To allow our horizons to evolve into non- 
symmetric ones, we would need to consider the general ver- 
sion of Eq. ([8]l |13|. 

With all of this mathematical background in mind, we turn 
to the Vaidya spacetime, which will be our chief example and 
testing ground for ideas. 



which we expect (and will see) to occur at large A: in this case 
"large" means large relative to nii. That is, setting 



A = miL 



(26) 



A is large if L ^ 1 . In later perturbative calculations we will 
use ( 1 /L) as the expansion parameter 

It is convenient to continue this transformation into dimen- 
sionless quantities. Defining a new radial coordinate R by 



III. VAIDYA BLACK HOLES AND THEIR HORIZONS 



Rm\ 



(27) 



In this section we identify and classify the horizons of the 
Vaidya spacetime. We start with a brief review of the space- 
time itself along with some assumptions specific to our calcu- 
lations. 



A. Vaidya spacetimes 

The Vaidya spacetime is described by the metric ll24l 

2m (v) ~ 



the metric becomes 



ds"- 



and so has stress-energy tensor 



dv^ + ldvdr + r^dd?- 



dm/dv 

Tab = —. — ij-[dv\a[dv\^ 



(21) 



(22) 



where m{v) is a mass function and dQ? = d9^ + sin^ 9d^^. 
Clearly if m(v) is constant, this reduces to the Schwarzschild 
spacetime in ingoing Eddington-Finkelstein coordinates and 
as for that case, v is an advanced time coordinate labelling 
families of ingoing radial null geodesies. Applying the energy 
conditions to ( |22] i, one finds that the mass function must be 
non-decreasing. In this case, the matter content of this space- 
time consists of infalling shells of null dust. 

For most of this paper we will be interested in Vaidya black 
holes that transition from an initial mass m; to a final mass 
m2. That is 



lim m{v)=mi and limm(v)=m2. 



(23) 



J 2 2 

ds = m-i 



-(l-^-^y^dV' + lLdVdR + R 



^dD} 



(28) 



It turns out that having the scaling parameter L appear directly 
in the metric in this way (rather than inside the mass function) 
is computationally convenient. 

For the numerical calculations we will need a concrete form 
for the mass function. In that case we almost always assume 
that 



M{V) = ^ + ^erf(y) 



(29) 



where erf is the error function. This describes a black hole 
irradiated with null dust that transitions from mass m\ to mass 
m2 — 2m 1. 

Next we locate and classify the horizons in Vaidya, begin- 
ning with the geometric horizons (which are easier). 



B. Geometric horizons 

Let us consider the geometry of the surfaces of constant v 
and r. First, a pair of (future-oriented) outward and inward 
pointing null normals are: 



^ = 



d_ 1 
_d_ 



1- 



2m(v) \ d 



and 



(30) 



We further assume that the evolution is characterized by a 
time-scale A so that 



m(v) = m\M ( — 
VA 



(24) 



where we have adapted m\ as a standard unit in order to write 
the mass function in a dimensionless form. In our study we 
will often also adopt dimensionless time 



V = 



(25) 



which lets us to study a range of spacetimes for a given mass 
function MiV), simply by changing the length scale A. In 
particular this allows for running to the near-equilibrium limit 



These have the usual rescaling freedom of £ ^ M and n — > 
n/h for any positive function h{v, r), but the given forms are a 
computationally convenient choice. 

By Q the corresponding expansions are 



(r — 2m) 
6(1,) = ^ ^ and 0, 



\n) 



(31) 



Thus two-spheres are trapped if r < 2m, untrapped if r > 2m 
and marginally trapped if r = 2m. The marginally trapped 



surfaces foliate a marginally trapped tube for which Eq. ( 16 1 

gives 



Jf„% 



1 



{t) 



4m{vf 



<0, 



(32) 




L = 200 




-20 



L= 1 




FIG. 2: Event and apparent horizons for three transitions from raj to mi = 2mi. For ease of visual comparison (particularly for L = 1 and 
1/100), the vertical axis is shown in terms of v/nii rather than V. As in FIG. [T] the apparent horizon expansion is caused by the infalling 
matter, here represented by the shading (with density proportional to shading intensity). For small L this also appears to be the case for event 
horizons. In contrast for large L, event horizon expansion is curtailed by the arrival of the matter shells. 



(since Gab£"n^ = 0) and as such it is a future outer trapping 
horizon. If dm/dv = it is null, isolated, and non-expanding 
but when dm/dv > it is spacelike, dynamical and expanding. 
These conclusions are confirmed by the evolution parameter 



which from Eq. (15 1 is 



C = 2- 



. dm 
dv 



(33) 



Several examples of these horizons are shown in FIG. [2] As 
expected they expand in response to infalling matter and do 
not expand in the absence of such stimulus. 

We can also check the slowly evolving conditions to con- 
firm that it is near-equilibrium for large L. For sufficiently 
large L we have 



C -0, 



2»" 



■^r^in) 



Gahrfri 



M \ 
L V m\M^ J 



1 
L 



M 
m\M^ 



4m]M^ 



< 



4-mfM^ 



and (34) 



|k'*'0(«)| 



where we have switched to dimensionless coordinates and a 
dot indicates a derivative with respect to V. In this regime the 
FOTH is slowly evolving. 



C. Event horizons 

Next we locate event horizons in Vaidya. The key is to re- 
call that the event horizon is necessarily a null surface and 
given the spherical symmetry of Vaidya, we expect this null 
surface to be similarly symmetric. From the metric, spheri- 
cally symmetric surfaces parameterized by (v(A),r(A), 0,0) 
are null if and only if they are tangent to one of the two null 
vector fields £ and n defined in ([30|. 

The surfaces tangent to n are clearly infalling with 
v=constant while the "outgoing" null surfaces are solutions 
of: 



dr 
dv 



= ; 1- 



2m(y) 



(35) 



Quotation marks are applied since for r < 2m (v) these "out- 
going" surfaces are also infalling (for example the curves la- 
belled A in FIG.[T]i. This is consistent with our earher obser- 
vation that for r < 2m{v), 9u'\ < 0. 

The event horizon must be an example of a surface tan- 
gent to i since those tangent to n all end up in the singuality. 
The trick is locating the correct null surface, and it is here 
that a knowledge of the future is important. Assuming that 
the black hole ultimately settles down (or at least asymptotes 
to) a Schwarzschild black hole of mass m2, the event hori- 
zon should similarly settle down at r = 2m2- Thus, one takes 
r(oo) = 2m2 as the (future) boundary condition. 

Generally this equation isn't solvable in exact form, how- 
ever it can be solved both numerically for the mass function 
( [29l ) and then perturbatively in the near equilibrium regime. 
We do this in the next two subsections. 



1. Numerical event horizons 



To solve ( 35 i numerically we must first deal with the future 
boundary condition: limy^oor(v) — 2m2. This can't be im- 
plemented exactly but it is sufficient to simply pick a large Vf, 
where m « 1112 and apply the boundary condition there. This 
works because, integrating backwards in time, null surfaces 
close to the event horizon exponentially converge on it. 

This convergence property can be seen in FIG.fllhowever it 
can also be demonstrated by a perturbative calculation (which 
in fact holds much more generally for perturbations about 
any given outgoing null geodesic). Consider a (spherically 
symmetric) null surface defined by r(A) = r£//(A)(l + ^ (A)) 
where r£//(A) is the true event horizon and |^(A)| ^ 1 per- 
turbs this to a nearby null surface. Then by ([35]), 
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where r„ — rEnivo) and (^„ = ^(vo). Now, making our usual 
assumption that the black hole started out in an equilibrium 
state with m = mi < 1112, this can be (crudely) bounded as 
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and so ^ — > (into the past) with a timescale « 4ot2/oti. As 
such, one doesn't need to get the final boundary condition ex- 
actly right: one just needs to pick a v for which m2 — m{v) is 
sufficiently small. In our examples v — lOA is more than suf- 
ficient: in that case OT2 — m(lOA) w 10^^^ and it asymptotes 
closer to the true event horizon with a "time" scale of Av =16. 
This kind of reverse integration is a standard way to find 
event horizons in numerical relativity |8 1. With the boundary 
condition taken care of in this way, it is then straightforward 
to solve (J35I1. Fourth-order Runge Kutta is sufficient and sev- 
eral examples of event horizons found in this way are shown 
in FIG. |2] For large scaling parameter L, the MTT and event 
horizon are close together (at least in the coordinate sense). 
Further they both appear to evolve in the same way, expand- 
ing only in response to infalUng matter. This is to be expected 
since this is a near-equilibrium regime. For small L the hori- 
zons are further apart and the presence of matter is seen to 
brake the expansion of the event horizon rather than drive it. 
Again this behaviour is to be expected, since in this case the 
full Raychaudhuri equation governs the evolution. 



2. Slowly evolving event horizons 



We now consider the large L limit. In terms of dimension- 



less quantities ( 35 1 becomes 
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This can be solved perturbatively in powers of \/L. To do 
this it is necessary to assume that over the range of interest, 
derivatives of M(y) are at most comparable in size to MiV) 
itself. 

For large A we expect the event horizon to be close to the 
geometric horizon and so consider solutions of the form 
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where we assume the coefficients (and their derivatives) are 
all much smaller than L. Then again representing derivatives 



of M with respect to V by dots, we can solve ( 38 1 to find 
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to second order. It is straightforward to extend this solution to 
higher orders, but for our purposes this is enough. 

At first glance there is something a bit strange about this 
solution: we started out solving for general "out-going" null 
geodesic solutions but in the end came up with just a single so- 
lution. The loss of other solutions can be attributed to the strict 



conditions imposed by the ansatz. As we saw in the last sec- 
tion, radial null geodesies rapidly converge to the event hori- 
zon in the past but correspondingly rapidly diverge into the 
future. Then it is our condition that derivatives of coefficients 
in the ansatz are comparable in size to the coefficients them- 
selves that eliminates all other candidates. We will return to 



this point in section IV C 2 where we perturbatively solve the 
general (radial) geodesic equation with a more general ansatz. 
That said, our sole solution is a good candidate for the event 
horizon: it is null to order 0{l/L^) and has the correct asymp- 
totic behaviour It may also be familiar to the reader as the 
event horizon candidate obtained in |25| (and also for linear 
mass functions in III6I ). In li25J the out-going null curve equa- 



tion ( 35 I was solved using a derivative expansion that is es- 



sentially equivalent to our procedure. There it was noted that 
this is an event horizon that is defined by the local geometry 
of spacetime. This is true but it is important to keep in mind 
that its identification as an event horizon depends on two as- 
sumptions: 1) at some point in the future m{v) -^ m2 and 2) 
the hierarchy of derivatives (or equivalently our ansatz) con- 
tinues to hold into the future. Thus claiming that it is a locally 
defined event horizon is slightly misleading. Event horizons 
are always null surfaces and the property of being null is lo- 
cal. The problem with locating an event horizon is picking the 
correct null surface and it is this that requires knowledge of 
the future. This case is no different. 

Finally we check the slowly evolving conditions ( 12 1. Ex- 



panding with (40 1 and using some results from section IIIB 
we find that to leading order in each quantity 
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on the perturbative event horizon. Thus for large A, 9h^ /2 <C 

Gahf"^'^ and \^ (9t£\ I <^ ^edm and this horizon is also slowly 
evolving. 



IV. PROPER TIME BETWEEN HORIZONS 

Finally we quantify the distance between the horizons and 
show that for large L the horizons become arbitrarily close. 
We begin with some general theory about hypersurfaces and 
their normal geodesies. 



A. Mathematical background 

A quick examination of FIG. [2] shows that for large L, the 
event and apparent horizons appear to be very close together 
while for small L they are relatively far apart. This has a good 
geometric interpretation: the r-values measure the areal radii 
of the horizons, so comparing these at constant v is actually 
comparing the areal radii of the horizons along the ingoing 
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K <0 

(geodesies converge) 
E. 



K >0 
geodesies diverge) 




FIG. 3: A spacelike surface £ in Minkowsld space along with some 
of its timelike geodesic normals and a conjugate point. The normal 
vectors are assumed to be past-oriented (as they are in our problem). 
By the theorem, the maximum proper time from qi to Z is measured 
along y. However, q2 is beyond the conjugate point x and so the 
maximum proper time from q2 to Z may or may not be along 7. 



radial null geodesies. This is a possible way to compare hori- 
zons at least in spherical symmetry and has recently been con- 
sidered in some detail in [.16.. 26 J . There it was shown that, for 
a range of matter models, in the near-equilibrium regime the 
two horizons will be close in areaf3T|. 

That said, this is not a measure of the distance between the 
horizons. Areal radius is a measure of area not distance and 
since curves of constant v are null, the geometric distance be- 
tween any two points along such a curve is zero. For a dis- 
tance we need a timelike or spacelike curve and in this case 
there is one obvious candidate: the timelike geodesic normal 
to the (spacelike) geometric horizon. With this candidate, the 
distance between the horizons is the proper time measured 
(backwards in time) along this geodesic from the geometric 
to the event horizon. The use of this proper time as a measure 
of distance is supported by the following theorem|T|: 

Theorem 1: Let ybe a smooth timelike curve connecting a 
point q (z M to a point p on a smooth spacelike hypersurface 
E. Then the necessary and sufficient condition that 7 locally 
maximize the proper time between q and E over smooth one- 
parameter variations is that y be a geodesic orthogonal to E 
with no conjugate point to E between E and q. 

Recall that given the congruence of timelike geodesies nor- 
mal to E with unit tangent vector field u", a point x on a partic- 
ular geodesic 7 is conjugate to E if there exists an orthogonal 
(Jacobi) vector field x" along 7 that satisfies the geodesic de- 
viation equation 



«"V„(«''Va")+^„/./«"z''«'' 



0, 



(42) 



and which vanishes at x but not on E. ^^ib/ ^^ ^^^ ^^^' 
mann tensor. Intuitively these Jacobi fields generate the vari- 
ations that map members of the congruence into each other, 
so a solution of these equations indicates the existence of two 
"nearby" members of the congruence that intersect at x. 

The theorem is demonstrated for the simple case of a 
Minkowski background in FIG. l3] which also demonstrates 
that it is a generalization of a Euclidean result: in (for exam- 
ple) three-dimensions the shortest distance between a point 
and a plane is measured along the straight line through the 
point that is perpendicular to the plane. The switch to timelike 



geodesies means that the shortest distance becomes a longest 
one, and allowing for curved E forces the inclusion of the 
qualification re conjugate points. 

As suggested by that diagram, the sign of the (trace of 
the) extrinsic curvature can imply the existence of conjugate 
points. This is the focusing theorem Q: 

Theorem 2: Let {M,gab) be a spacetime satisfying 
l%ah^"^ > Q for all timelike £," and let T, be a spacelike hy- 
persurface with K = Ou\ < at a point p £'L. Then within 
proper time T < 3/\K\ there exists a point q conjugate to E 
along the geodesic 7 orthogonal to E and passing through p, 
assuming that 7 can be extended that far 

Paraphrasing, this theorem says that initially converging 
geodesies K < necessarily converge in finite time. In fact 
it can be seen quite easily. For a hypersurface orthogonal con- 
gruence of timelike geodesies u", the timelike Raychaudhuri 
equation tells us that 
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(43) 



In our case the geodesies start out being orthogonal to the 
spacelike horizon which has induced metric qij = efe'^-g^b and 

extrinsic curvature Kjj = e"e^V„Mi (e" is the pull-back oper- 
ator onto the surface). Then at the point that they cross the 
horizon, the geodesies have expansion 0(„) — K — q'-'Kij and 

shear ff; ■ — Kij — \Kqij. Thus if the dominant energy condi- 
tion holds {^abu"u'' > 0), it is trivial that 
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dz - 3 
and the result follows by direct integration. 

B. Finding the normal geodesies 

In order to apply the preceding theory, we first need to find 
the normal geodesies and it is to this problem that we now 
turn. 

Parameterizing with the non-affine (but otherwise very con- 
venient) parameter v, the radial geodesic equations for the 
Vaidya spacetime are: 
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where X"{v) = [v, r{v) ,0,0] and on the right-hand side the ex- 
tra function / is required since v is non-affine. From the first 
equation (X" = v) we immediately find that / — FJ^^, and so 
general radial geodesies are solutions of 
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(47) 



are the relevant non-vanishing Christoffel symbols. This is 
a general equation for radial geodesies and so solutions can 
be timelike, spacelike or null. In particular the event horizon 
will be a null solution to these equations (in this case they 
are equivalent to Eq. 35 1. The parameterization by v means 



that certain other solutions will be excluded: most notably the 
ingoing v ^constant null geodesies. However, such solutions 
are not of interest to us in what follows. 

The timelike normal geodesies are solutions of these equa- 
tions. Their initial conditions are determined by the time- 
like normal to the spacelike geometric horizon. The future- 
oriented version of this normal is 
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where here m = dm/dv. Then the initial conditions for the 
timelike normal geodesies are 
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These equations are not solvable in closed form however 
it is straightforward to solve them numerically. We will be 



mainly interested in the mass function (29 1 and in particu 



lar studying how geodesic properties change with the scale 
parameter L. Then, as for the null geodesic equation, it is 
numerically and theoretically advantageous to switch to di- 
mensionless form so that the mass function is invariant with 
respect to L and the geodesic equations become 
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while the initial conditions maintain their original form 
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FIG. 4: The geometric horizon and its timelike geodesic normals for 
L = 1 . The lower part of the figure is an enlargement of the boxed 
region in the upper part of the figure. Nearby geodesies intersect out- 
side the event horizon, suggesting the existence of conjugate points 
there. The dashed part of the geometric horizon has A"/ „) < (past 
oriented-normal) . 



For L < 100, a standard technique such as fourth-order 
Runge-Kutta is sufficient to solve the equations over much 
of the range of interest. However for larger values of L or 
|y I > 5 there are numerical difficulties. In both of those sit- 
uations we enter a regime where all of the curves of inter- 
est are closer together than any reasonable numerical accu- 
racy that we might chose to work with. In such cases we 
turn to perturbative methods and consider the regime where 
R = 2M(y)(l +p(y)) with p{V) < 1. The details of how 
that expansion is done differs for the two cases and so we 
leave the details of this expansion to the appropriate sections. 



C. Spacetimes with L < 100 

1. Numerical regime 

We begin with the regime where the geodesic equations can 
be solved numerically. In particular we would like to know 
when (if) the normal geodesies intersect the horizon and what 



is their proper time length; as set out in section pV A| this will 
give a measure of distance between the horizons. 



To study these issues, we numerically solve equation ( 35 1 
with the appropriate (future) boundary condition to locate the 
event horizon. We then step along the geometric horizon re- 



peatedly solving (50 1 with boundary conditions (51 1 to find 



the normal timelike geodesies. Next, we search for intersec- 
tions between the geodesies and event horizon. The Newton- 
Raphson method turns out to be sufficient for this purpose; 
if the numerics are reliable and there is actually is an in- 
tersection, experience shows that this method quickly con- 
verges. With an intersection in hand, proper time length is 
easily found by integrating 
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where Vf is the value of v at which the geodesic departs 
the geometric horizon and v,„, is the value where it inter- 
sects the event horizon. If there are no conjugate points 
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along these geodesies, then by Theorem 1 the time along any 
given geodesic will (at least locally) be the maximum possi- 
ble proper time that can be measured along any timelike curve 
connecting the intersection point to geometric horizon. 

This process is demonstrated in FIG. [4] for the usual error- 
function-determined mass with L = I. The top half of the 
figure shows the normal geodesies over a wide v range while 
the bottom half magnifies the region over which the geometric 
horizon is dynamically evolving. Keep in mind that this is 
a general coordinate system and so one should not jump to 
any conclusions based on the apparent shape of curves. For 
example in this coordinate system the normal geodesies do not 
look like they are orthogonal to the geometric horizon (though 
of course they are) and the extrinsic curvature of the geometric 
horizon is not obvious. For future pointing normal u", the 
trace of the extrinsic curvature is 



%)-^"''V„Mfc = 



m — 2min — 8m 
8mm^/^ 



(53) 



where in this case with L=l, dots are regular derivatives with 
respect to v. The sign of this quantity is shown in the bot- 
tom half of the diagram and provides a nice demonstration 
of the focussing theorem. By that theorem, conjugate points 
should exist along normal geodesies which exit the geometric 
horizon when ^(_a) = — ^(«) < (the negative sign is intro- 
duced since we are interested in past-oriented geodesies). A 
line of of these conjugate points can clearly be seen in the 
bottom half of the figure and indeed these do occur on the 
expected geodesies, though outside of the event horizon. For 
the region over which ^(-„) > 0, where the theorem makes no 
predictions, there do not appear to be any conjugate points at 
all. Thus in this example, independent of the curvature of the 
apparent horizon, the normal geodesies maximize the proper 
time between the geometric and event horizons. 

All other values of L that we examined also lacked conju- 
gate points between the horizons and so the geodesies in those 
cases continue to maximize proper time between the horizons 



(though see section IV E for a more complicated situation). 
Three representative examples are shown in FIG. |5] Other 
trends apparent in those figures also continue for other values 
of L: for L < 1/10 the apparent horizon transition becomes 
sharper but otherwise the figure is more-or-less unchanged 
while for L > 10 the curves become closer together on the fig- 
ure and the positive/negative v behaviour becomes more sym- 
metrical. 

It FIG. PI another important trend is clear: as Vf decreases, 
the geodesies increasingly asymptote towards the event hori- 
zon rather than clearly intersecting it. This effect can also be 
seen in FlG.fflwhere for larger values of Vf the geodesies cross 
the event horizon and then appear to be asymptoting "back" 
towards it from the outside, while for more negative values 
they instead appear to be asymptoting from the inside and it is 
not at all obvious that they ever intersect. The non-existence 
of an intersection in such a case is, of course, very difficult 
to show numerically: however far we track the curves back- 
wards in time (and so to whatever numerical accuracy) it is 
always possible that an intersection occurs just a bit further 
back (and so beyond any accuracy that we choose). Thus to 



demonstrate the existence of such non-intersecting geodesies 
we turn to perturbative calculations. 



2. Perturbative regime 

For negative V, we can asymptotically expand the mass 
function (|29ll and its derivatives as 
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Then, for example, ify = — 10, M— 1 and its derivatives are 
all less than 10^'*''. Meanwhile, as even a cursory examination 
of FIG. |4] and FIG. |5] demonstrates, if Rjq is the location of 
the timelike geodesic then Rtg ^ 2M will generally be much 
larger In such a regime it is reasonable to take M = 1 and set 
its derivatives to zero; essentially we revert to Schwarzschild 
geodesies. 

Then we look for perturbative solutions, expanding around 
/? = 2(thatisr = 2mi): 

Rtg = 2{\+p[V)). (55) 

To first order in p the geodesic equation ([50| becomes 



1 cP-p 3 dp p 
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which has general solution 
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for some constants A^ and Bo. Going back to the metric, it is 
easy to see that the value of B^ determines the signature of the 
geodesic: Bg > 0^ spacelike, £„ = <=> null, and B„ < <^ 
timelike. 

Next, in many of the numerically ambiguous cases, the 
regimes where the numerical and perturbative analyses work 
overlap and so we can use these perturbative solutions to study 
the ultimate fate of the geodesies. For a Vo in this overlap (typ- 
ically — 10 < y„ < —40 is used in our work) we can determine 
the coefficients in the perturbative solutions by matching to 
the numerics. Doing this 
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p' ^ =(2p„' ^ - 4p„' ^/L) exp - (y - Vo) (59) 
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where Po — p{Vo) and Po = dp/dV\ 
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FIG. 5: The timelike geodesies normal to geometric horizons for three representative values of the scale parameter L. Again the shading shows 
the concentration of inf ailing matter and the vertical axis is v/mj rather than V. To reduce clutter, the geodesies are cut off at the event horizon. 



Now, if no intersection has occurred for V > Vo, we can use 
these matched approximate solutions to study potential V < V„ 



intersections. First note that using (58 i we can rewrite (59 1 as 
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By the discussion following ( 57 1, p ™ — 4pJ'^ > and so we 
see that if a timelike geodesic arrives at Vo with 
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(61) 



and our approximation holds, then it will never intersect the 
event horizon. Numerical checks show that for L—\ this con- 
dition is satisfied for all geodesies which intersect the geomet- 
ric horizon at Vj < —1.087 and not satisfied for larger values 
of Vf. Similar results are found for the other values of L and 
so this suggests the following picture: in each case there will 
be a transition point which separates normal geodesies which 
intersect the event horizon from those that don't. 

In this situation one might guess that in the approach to 
the transition point from above, the intersection point diverges 
with Vint -^ ~°°- The perturbative solutions support this hy- 

gives 



pothesis. Solving (|59|) for p^^ = p^" 
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The denominator is always positive and the numerator van- 



ishes exactly where the bound in condition (61 1 is saturated. 
Thus, as expected, Vint -^ ~°° as Vf approaches the transition 
point. 



3. Proper time lengths 

With this picture in mind we turn back to the numerics 
and results that were obtained for intersections and geodesic 



il 

ill 

il 1/ \/ \ 






\2 







T 

mi 



\ 



FIG. 6: Proper time lengths of normal geodesies versus departure 
point from the geometric horizon for seven values of the scale pa- 
rameter L. Starting from the left-most and sharpest peak and ending 
with the right-most and flattest peak they are L = 10^ , L = 10^ , 
L = \0-\ L = 1, L = IQi and L = 10^, L = 10^ For V > Vpeak 
geodesies intersect the event horizon and their length between the 
horizons is show. For V < V„p„j. there is no intersection with the 
event horizon. Instead the total length of the geodesic from V = — oo 
to Vf is plotted. 



lengths. These results are summarized in FIG. [6] which re- 
quires a little explanation. First, for L < 10 there are sharp 
peaks with discontinuous first derivatives (there is also one for 
L — 100, though it cannot be seen well in this figure). Those 
peaks are the demarcation between geodesies that ultimately 
intersect the event horizon (on the right) and those that don't 
(on the left). Computationally, they correspond to the point 
where the intersection finder first fails. However, within nu- 
merical error they also match the transition points predicted 
above. To the right of the peak the length from V,,,, to Vf is 
plotted however on the left it is instead the total length of the 
geodesic from V = —°° to Vf (see Appendix[A|for a comment 
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on how this is done), hituitively one can think of these lengths 
as arising from a competition between how close a curve is to 
being null versus its coordinate length. Thus for Vf < Vpeuk all 
curves have infinite coordinate length but are much closer to 
being null as Vf decreases (have another look at the geodesies 
in FIG. ph. Meanwhile as Vf increases past Vpeak the coordi- 
nate length of curves decreases and they ultimately become 
"more null" and so their decrease in length is not unexpected. 

From the point of view of defining a maximum distance be- 
tween the event and geometric horizons it is the T {Vpeak) that 
is important. Assuming that there are never any conjugate 
points between the horizons (we have checked many specific 
cases but not found a general proof), then Theorem 1 tells 
us that from a given point on the event horizon, the (at least 
locally) maximum proper time from that point to geometric 
horizon is along a normal geodesic to the geometric horizon. 
This is exactly what is plotted in FIG. l6] though by the inter- 
section with the geometric, rather than event horizon. Since 
the peak corresponds to V = — oo on the event horizon, it is 
only the part of the curve Vf > Vpeak that is relevant to this dis- 
cussion. Then T (Vpeak) is the maximum proper time distance 
between the horizons and this distance is measured along an 
ingoing timelike geodesic that asymptotes to the event horizon 
at y = — oo and ultimately orthogonally intersects the geomet- 
ric horizon at Vpeak- 

The second important observation to make from FIG. l6]is 
that as L increases, that maximum proper time decreases. This 
is slightly qualified by the fact that numerical difficulties arise 
for L = 1000 as the coordinate values of the curves become 
extremely close together; in fact for L = 1000 the algorithm 
doesn't find any intersections and so all integrals are estimated 
using the methods of Appendix [A] That said, the trend appears 
clear and this gives support to the hypothesis that for large 
enough L the horizons become arbitrarily close. We will close 
the gap in this argument in the next section where we will 
employ perturbation techniques to study this large L limit. 



coefficient of the differential equation) as: 



D. Large L spacetimes 

We now consider a regime where the entire geodesies can 
be described perturbatively. In particular we confirm the con- 
jecture that in the slowly evolving limit (as L — > oo) the maxi- 
mum distance between the geometric and event horizons goes 
to zero. 



1. Geodesies 
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where M = dM/dV and M = d^M/dV'^. 

To obtain a general solution of the geodesic equation we 
have to allow derivatives of p to be much larger than p itself. 
We have already seen this kind of behaviour during our first 
venture into perturbative solutions ( [57] i. Inspired by that work 
but realizing that M{V) is no longer constant, we begin by 
trying a first order solution of the form 
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where /i/4(V) and /i/2(V) are functions of V and all of the L 
dependence is contained in the 
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(this reduces to (57i for M ~ 1). For normal geodesies the 



"initial" value in the integral will usually be the departure 
point from geometric horizon and so it is written as Vf. Substi- 
tuting, it is straightforward to see that (64 1 in its general form 



is a first order solution to ( 63 1. This can be understood by not- 
ing that to lowest order, the /i m and /1/2 will be effectively 
constant relative to X (for which changes in V are magnified 
by the L). 

This freedom in the first order solution is eliminated by go- 
ing to second order in the differential equation. Then, noticing 
the p2 terms in (63 1, we try a solution of the form: 



P2=P\ + j2 (32m2 + 16MM) 



(66) 



This is a second order solution if and only if: 

■^1/4 - Jp ' -^1/2 - ^ ' ^3/4 - --f^ (67) 

92 



and gi 



K 

M8' 



To this end we return to the geodesic equation (50i. This for arbitrary constants A^ and B„ and functions giuiV) and 
time we assume p ^ j^<^l and expand to third order (in each giii^Y)- 
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Just as we had to go to second order to finalize the first 
order solution, we go to third order to finalize the second order 
solution. Then 



h/4 ■ 



h/2 ■ 



1 

1 
'M4 



a„ 



-AA, 



(2MM + M^- 



"kV 






M 
-2B„ 



dv 



(68) 



3M- 



2M2 
If 



where a„ and j3o are new constants and we omit the informa- 
tion on third order quantities since they won't be required in 
what follows. 

This expression for general geodesies is quite formidable 
however we can extract useful information from it. First note 
that 



exp 



4(m2/mi) 



<X < exp 



4 



(69) 



and so with L large, all of these solutions will rapidly asymp- 
tote to Peh for V < Vf. Further they will rapidly diverge (and 
leave the regime of validity of the approximation) for V >Vf. 
This helps to explain why Peh contains no free constants: it 
is the only solution that we have found that stays close to the 
geometric horizon. Other geodesies, whether timelike or null, 
diverge on a time scale of 1 /L. 

Now we focus our attention on the timelike normals. The 
initial conditions ( |5T] i become 



P{Vf)^0 and 



dp 
dV 



Vf 



2% 
Mf 



(70) 



where Mf — M{Vf) and Mf = M{Vf). Applying these we find 
that 



Aa = and B„ = -AMfMJ , 



(71) 



while 



finding where the first Une of (73 i vanishes. That is we need 
to solve: 



■ {MfMf 



+ 2M^ 



MfMf 



1 (r ^ dy 

M{V) ^^^ \Jvf 4M(v) 



(74) 



If M{V) is positive, bounded from below and increasing 
(which it is) then it is clear that this can only have solutions if 
Mf is negative or (as a stronger condition) 



Mf< 



2MJ 



M 



(75) 



/ 



Note too that if this bound holds, then there will always be a 
solution since the right-hand side of equation (|74]i is monoton- 
ically increasing and goes to zero as V ^- — oo. Equivalently, 
when the bound is saturated we will necessarily find the inter- 
section value Vint -^ — °o. 

Based on our experience in the previous section we ex- 
pect such a transition to correspond to a peak in the normal- 
geodesic proper-time length functions. Indeed this seems to 
be the case. For our particular mass function the bound is 
saturated when 



2 "^ 

(3+eify)y--'^=o, 



(76) 



which has solution V w 0.30765. Examining FIG. [6] shows 
that indeed for L = 100 and L — 1000, the peak in the length 
function does occur at about this point. 

So the pattern seen in the last section continues into this 
regime: for our mass function, there is a Vpeak which divides 
the normal geodesies that intersect the event horizon from 
those that never intersect it. Let us now consider the proper 
time lengths of those geodesies. It will be sufficient to work 



at lowest order. First, from the metric ( 28 1 the (square of the) 



rate of change of proper time along these geodesies is 



Uo = -16MJ{mfMf-2MJ) and j3„ = 16MJM] 



IflVlf . 



(72) 



Then, the timelike normal geodesic that intersects the geomet- 
ric horizon at Vf is parameterized (to second order in 1 /L) as 



Ptg —Peh 



AMfMJ 

M4 



X' 



16M: 



/MfMf + 2MJ\ X 




l6MfMJ 

/ IGMJMf ' 
I M8 



Mf 



3M + 2- 



.M^ 

1a 



dv 



L2 






This is still not a simple expression, however, again it is 
useful. From our experience with the numerical results, we 
expect any intersections with the event horizon to occur when 
[Vf — Vi„,)L ^ 10. If this is true then X,>„ will be small and so 
we should be able to obtain a good approximation of V,-,,, by 



-f-V 

ffij \dv J 
So, substituting in 



'PtgL'-~AL[M 



dp 



TG 



dV 



Ptg ' 

we find 

1 d% 



AM AMf (M'f 

77 "^^tT I m4 



X' 



(77) 



(78) 




mxdY ' V f \m'^ ] ''"*' V^v;;- 4M(v) 



Now the proper length from the geometric horizon to the event 
horizon will be bounded by the integral along the geodesic 
from Vy to — oo 



AJ I-'" ( dx\ 
mi Jvf \dV J 



(80) 
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and in turn we can use the fact that M > 1 to bound this inte- 
gral by 



(81) 



Thus, for any non-decreasing positive mass function 
m(v) = miM{\'/L), the proper time distance between the hori- 
zons goes to zero as L — ?► oo. In this case the event and geo- 
metric horizons are geometrically close. 



E. More complex spacetimes 

As a final example we consider a more complicated mass 
function in which two distinct pulses of dust fall into a pre- 
existing black hole: 



m{v) 



mi(3+erf(v)+erf(v-20)) 
2 



(82) 



The resulting spacetime along with its horizons and geodesies 
is shown in FIG. |7] 

Turning to the proper time graph, this time there are three 
peaks and each of these corresponds to a transition between 
geodesies intersecting the event horizon and geodesies which 
diverge to -~°°. The first and third peaks correspond to those 
seen in the earlier examples, while the middle peak arises as 
some normal geodesies sneak through the gap where the hori- 
zons almost meet and then head off to v = — oo. While the 
proper time curve continues to record local maximum dis- 
tances, it is clear that this time it does not always record global 
maxima: the same points of intersection with the event hori- 
zon are necessarily recorded repeatedly in the lead-up to each 
of the three peaks. 

Finally note that if we rescaled this mass function as 



m(v) = 



mi(3+erf(^)+erf(f-20)) 



(83) 



and considered the large L limit, all of our near-equilbrium 
results from the earlier sections would continue to hold. 



V. DISCUSSION 

From an astrophysical perspective one can reasonably ar- 
gue that attempting to define the exact boundary of a black 
hole is physically about as meaningful as discussing the num- 
ber of angels that can dance on the head of a pin. If one is 
only interested in what can be observed from afar then it is 
the near-horizon fields of the membrane paradigm that are the 
significant quantities. By construction the horizons are all un- 
observable. While this is certainly true in principle, in prac- 
tical calculations this apparent irrelevance often evaporates: 
both numerically and perturbatively horizons of all types are 
often observed reacting to their environment in intuitive ways. 
The reason for this is that while some particularly dramatic 




FIG. 7: Two shells of dust fall into a black hole causing two distinct 
phases of expansion. As in the previous figure the top graph shows 
horizons and timelike normal geodesies while the bottom graph is of 
the proper time between horizons measured along those geodesies. 



events (such as the actual merger phase of a black hole col- 
lision) are far from equilibrium, most black hole physics is 
near-equilibrium. In fact as discussed in l;15j even some very 
dramatic events, such as black hole formation, can be in this 
regime. Then it is reasonable to expect the horizons to be 
"close" together, jointly reflecting the evolution of the exter- 
nal geometry. 

In this paper we have begun to quantify this idea. For the 
Vaidya spacetime we have located event and apparent hori- 
zons and measured the distance between them using the time- 
like geodesies normal to apparent horizon. We have seen 
that at most these proper-times are of the order of the mass 
of the black hole while in the near-equilbrium regime time 
separation goes to zero. If the black hole is eternally near- 
equilibrium the horizons are close together. If there are peri- 
odic bursts of far-from-equilibrium activity, the horizons will 
separate. However, once things begin to calm down they come 
together again. 

Of course, all of this was done for spherically symmetric 
surfaces in spherically symmetric spacetimes and those are 
not real world conditions: many nice ideas in general relativ- 
ity have foundered during their generalization from spherical 
symmetry. That said we believe that a generalization in this 
case is quite feasible. It should be possible to perturbatively 
construct the spacetime geometry around a slowly evolving 
(apparent) horizon and we expect that in that geometry it will 
always be possible to locate an event horizon candidate: a null 
surface that hugs the apparent horizon. Such surfaces have al- 
ready been located and studied in studies of black brane space- 
times Il25ll27ll28ll . More speculatively one might also hope 
that this surface will (partially) resolve the non-uniqueness for 
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apparent horizons. It may turn out that there is a unique null 
surface close to all near-equilibrium apparent horizons. An 
investigation of these ideas will be the subject of a future pa- 
per 

There is one obvious caveat in our arguments that the ex- 
act location of the boundary of a black hole may not matter. 
For stationary black holes it is well-known that the area of 
the event horizon is directly proportional to the entropy. Fur- 
ther, we have seen that the second-laws of horizon mechan- 
ics show that horizon area, like entropy, is non-decreasing in 
time. Thus, it is commonly assumed that entropy continues 
to be proportional to horizon area for dynamical black holes. 
Then it does become important to pick a single horizon so that 
entropy will be uniquely defined (see lfT6ll for further discus- 
sion on this point). 

Note however that there is another possibility. It is well- 
known that there is no local definition of gravitational energy 
in general relativity. There may, or may not, be a unique defi- 
nition of quasilocal energy 1 29 1 in stationary spacetimes, how- 
ever far from equilibrium its definition is even more problem- 
atic. Now, thermodynamic entropy is operationally defined 
by the entropy flow equation 5E = T5S. Thus if energy be- 
comes ill-defined or non-unique (not to mention temperature) 
it would not be surprising if there were similar difficulties for 
entropy. As long as any definition of entropy settles down to 
the correct value in equilibrium this may be enough: since all 
the horizons join (or at least asymptote) together in equilib- 
rium any one of them would be a reasonable candidate. Many 
standard notions of physics (including of course the flow of 
time) are drastically altered in strong gravitational fields and 
it is quite possible that entropy might join the list - partic- 
ularly in the regime of non-equlibrium thermodynamics. As 
long as there are no violations of our standard laws of physics 
for observers in the weak-field regime, complications in the 
strong-field regime are acceptable. Some further discussion 
of these ideas may be found in 1281 which tries to reconcile 
notions horizon-defined entropy with those arising from the 
AdS-CFT correspondence. 
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Appendix A: Measuring infinite geodesies 



A final application of ( 59 1 arises when we wish to find a 
bound on the total length of a geodesic normal between ¥„ 
and V = — oo. This situation arises both in cases where the 
geodesic never intersects the event horizon but instead contin- 
ues on between the horizons to V = — oo and also in the case 
where it may intersect for some value of V < Vo but the nu- 
merical precision isn't high enough to find that intersection. 
In the second case the proper time from V = —°° to V„ will 
still provide an upper bound on the true length. 

We proceed as follows. By the metric (28 1 the rate of 



change of proper time along a timelike curve in the asymp- 
totic regime (where p is small and derivatives of m can be 
neglected) is 



1 fdTV 
dv) 



^ dV 



(Al) 



and so by Eq. ( 59 1 



dx 
dV 



' m\L 



\ 



pTG_ 



iTG 



dV 



expl ^(V-K)) .(A2) 



Then the proper-time length of the timelike geodesic for V < 
Vu (whether or not it intersects the horizon) is bounded by 



Ar = 4m 



pr-^Pi' 



(A3) 



The values for pj'^ and pj*^ can be calculated from the 
numerical simulations. For the examples considered in FIG. [6] 
we find Ar < lO^^mi is negligible compared to the proper 
time along the geodesic from Vo to Vf. As such in FIG. [6] 
proper times to the left of the peak are estimated by integrating 
from Vo to Vf. 



[1] P. Anninos et. al, Phys. Rev. Lett. 71 2851 (1993); P An- 

ninos et. al., Phys. Rev. Lett. 74, 630 (1995); Baiotti et.al, 

Phys. Rev. D 71 024035 (2005); Schnetter et. al. Phys. Rev. D 

74 024028 (2006) 
[2] E. Poisson, Phys. Rev. D 70 084044 (2004); E. Poisson, 

Phys. Rev. Lett. 94 161103 (2005); E. Poisson and I. Vlasov, 

Phys. Rev. D 81 024029 (2010). 
[3] K. S. Thome, R. H. Price and D. A. Macdonald, Black holes: 

the membrane paradigm, (Yale University Press, 1986). 
[4] L. RezzoUa, R. P. Macedo and J. L. Jaramillo, Phys. Rev. Lett. 

104,221101 (2010). 
[5] R. Penrose, Phys. Rev. Lett. 14 57 (1965). 



[6] S.W. Hawking and G.F.R. Ellis, The large scale structure of 
space-time, (Cambridge University Press, 1973). 

[7] R. M. Wald, General Relativity, (University Of Chicago Press, 
1984). 

[8] J. Thornburg, Living Rev. Rel. 10, (2007) URL: 
http://www.livingreviews.org/lrr-2007-3 (cited on Septem- 
ber 10, 2009) 

[9] S.A. Hay ward, Phys. Rev. D49 6467 (1994). 
[10] A. Ashtekar, S. Fairhurst and B. Krishnan, Phys. Rev. D62 
104025 (2000); A. Ashtekar, C. Beetle and J. Lewandowski, 
Phys.Rev. D64 044016 (2001); A. Ashtekar, C. Beetle, and 
J. Lewandowski, Class. Quant. Grav. 19 1195 (2002); 



16 



[11] A. Ashtekar and B. Kiishnan, Phys. Rev. Lett. 89 261101 
(2002); Phys. Rev. D68 104030 (2003). 

[12] A. Ashtekar and G. Galloway, Adv. Theor. Math. Phys. 9 1 
(2005). 

[13] I. Booth and S. Fairhurst, Phys.Rev.D 75:084019 (2007). 

[14] E. Schnetter and B. Krishnan, Phys. Rev. D 73 021502 (2006). 

[15] W. Kavanagh and I. Booth, Phys. Rev. D 74 044027 (2006). 

[16] A. Nielsen, The spatial relation between the event horizon and 
trapping horizon, arXiv: 1006.2448 (2010). 

[17] S. W. Hawking and J. B. Hartle, Commun. Math. Phys. 27 283 
(1972). 

[18] E. Gourgoulhon and J. L. Jaramillo, Phys. Rev. D 74:087502 
(2006). 

[19] A. Ashtekar and B. Krisnan, LivingRev.Rel.7:10 (2004) URL: 
http://www.livingreviews.org/lrr-2004-10 (cited on February 
17, 2010); I. Booth, Can.J.Phys. 83 1073 (2005); E. Gourgoul- 
hon and J. L. Jaramillo, Phys. Rep.423, 159 (2006). 

[20] I. Booth, L. Brits, J. Gonzalez, and C. Van Den Broeck, 
Class. Quant. Grav. 23, 413 (2006). 

[21] J. L. Jaramillo et. al, AIP Conf. Proc. 1122, 308 (2009). 

[22] T. Chu, H. P. Pfeiffer and M. I. Cohen, Horizon dynamics in 
perturbed Kerr spacetimes, arXiv: 101 1.2601 (2010). 

[23] I. Booth and S. Fairhurst, Phys. Rev. Lett. 92 01 1 102 (2004). 



[24] E. Poisson, A Relativists Toolkit: The Mathematics of Black- 
Hole Mechanics, (Cambridge University Press, 2004). 

[25] S. Bhattacharyya et al, JHEP 0806, 055 (2008). 

[26] A. B. Nielsen, M. Jasiulek, B. Krishnan and E. Schnetter, The 
slicing dependence of non-spherically symmetric quasi-local 
horizons in Vaidya Spacetimes, arXiv: 1007.2990 (2010). 

[27] S. Bhattacharyya, V. E. Hubeny, S. Minwalla and M. Ranga- 
mani, JHEP 0802, 045 (2008). 

[28] I. Booth, M. Heller, M. Spalinski, Phys. Rev. D 80, 126013 
(2009). 

[29] L. Szabados, Living Rev. Rel. 7 4 (2004), URL 
http://relativity.livingreviews.org/Articles/lrr-2004-4. 

[30] Alternatively, and perhaps more logically, this should be inter- 
preted as a Clausius-like statement of the second law: 5£t\fq ~ 
{VQGab(-''(-'')/Kf in analogy with SS > SQ/T. See 1 18| a dis- 
cussion of this point of view. 

[31] There is a difference in emphasis between 1161 and our paper. 
Here we are interested in the fractional difference between the 
horizons while |16| is more interested in the absolute differ- 
ence and the implications that difference might have for entropy 
and thermodynamics in general. In absolute terms the differ- 
ence may be quite large. 



